Relativity Restored: Dirac Anisotropy in QED 3 
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We show that at long lengthscales and low energies and to leading order in 1/N expansion, the 
anisotropic QED in 2+1 dimensions renormalizes to an isotropic limit. Consequently, the (Eu- 
clidean) relativistic invariance of the theory is spontaneously restored at the isotropic critical point, 
characterized by the anomalous dimension exponent of the Dirac fermion propagator n. We find 
n = 16/37r 2 iV. 



Quantum electrodynamics in (2+1) dimensions 
(QED3) has recently emerged as a low-energy effec- 
tive theory of a number of condensed matter systems 
0, §, |, |, §, |, 0. Examples range from fluctuating d- 
wave superconductors in uderdoped cuprates Q, 0, 0] to 
pyrolitic graphite @] to Hcisenberg antiferromagnets and 
spin liquids [0, 0] . While these multiple reincarnations 
of QED3 differ mightily in their physical content, they 
all share certain important formal similarities. The low 
energy behavior is controlled by an infra-red fixed point 
where the gauge field acquires a universal dimensionless 
coupling constant g oc 1/N, N being the number of 
Dirac fermion flavors || 0|. At values of g larger than 
some critical value g c (N < N c ) it is believed that the 
theory has an instability into a state with broken chiral 
symmetry, with fermions spontaneously acquiring a 
finite dynamical mass ||. Among the formal aspects 
shared by the above theories surely one of the most 
ubiquitous is the spacetime anisotropy - such low-energy 
effective theories are obviously only pretending to be 
"relativistic" . They hail from non-relativistic quantum 
Hamiltonians and are not obliged to be invariant under 
Lorentz transformations. Consequently, they often 
contain more than one "speed of light" resulting in the 
above anisotropy. An important question is to what 
extent are the properties of these effective theories sim- 
ilar to the genuine, isotropic QED3 and, in particular, 
what is the nature of the critical behavior and chiral 
symmetry breaking when such anisotropy is present. 

In this Letter we address the problem of Dirac 
anisotropy in QED3. Our point of departure is the as- 
sumption that the symmetric (massless or critical) phase 
of isotropic QED3, obtained when the number of fermion 
flavors is larger than a critical value, N > N c , is con- 
trolled by a stable non-trivial infrared critical point [0, 0] 
characterized by the anomalous dimension exponent r\ > 
of the gauge-invariant Dirac fermion propagator. Based 
on this assumption, which is almost certainly correct for 
N ^> N c , we derive the following results within a 1/N 
expansion: i ) When anisotropy is turned on at this inter- 



acting critical point we find it to be marginally irrelevant 
in a perturbative sense. This implies that the symmet- 
ric, critical phase of QED3 remains unaffected by small 
Dirac anisotropy. In particular, the value of r\ remains 
unchanged, ii) Going beyond the perturbative regime, 
and by exploring the structure of renormalization group 
(RG) flows, we argue that any finite anisotropy is also 
irrelevant. Finally, Hi) we compute the explicit value of 
n and find v = 16/3ir 2 N. Our results imply that the rela- 
tivistic invariance of a QED 3 -like effective theory is itself 
an emergent property: it is spontaneously dynamically 
restored at the critical point. 

The anisotropic QED3 can be defined as follows: 
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C = + ia„)W n) + ^(9x a) 2 (1) 



where ip^ n ' is a Fermi field associated with a node n, 7^ 
is a Dirac matrix, and a M is a massless U(l) gauge field 
related to fluctuations of unbound 2 + 1 vortex loops 



We also introduced the diagonal "nodal" metric g^j 
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forms can be reduced to this one by suitable rescalings 
of spacetime coordinates and fermion and gauge fields. 

The Dirac anisotropy of C (0) is more sinister that its 
bosonic kin JToj] . In the Higgs-Abelian gauge theory the 
anisotropy can be fully rescaled out of the matter part 
of the action leading to the new effective action with 
the anisotropy stored only in the gauge field Maxwellian 
action. Since matter cannot generate any anisotropic 
contribution to the gauge field and by the virtue of the 
isotropic charge being a relevant operator, it is easy to 
show that anisotropy of the bosonic theory is maginally 
irrelevant. 

In the fermionic QED3 (0) the above simple procedure 
does not work because one cannot simultaneously rescale 
the kinetic energy for all fermion species. We therefore 
keep the anisotropy confined to the matter part of (0) and 
proceed from there. The two-point vertex function of the 
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non-interacting theory for, say, (1,1) Dirac fermions is 

p(2)/ree = + + w A72 fc 2 (2) 

and the corresponding Green function equals 



GlS(k) 



(3) 



In what follows we assume that both vp and va are di- 
mensionless and that eventually one of them can be cho- 
sen to be unity by the appropriate choice for the "speed 
of light". The anisotropy parameter up, — vf/va =/= 1 
breaks the Lorentz invariance of the theory ([!]). However, 
the theory still respects time-reversal and parity and for 
-/V > N c the system is in its chirally symmetric phase 
jn]]. These symmetries force the fermion self-energy of 
the interacting theory to assume the following form: 



(4) 



The coefficients Q are in general different from unity. 
Furthermore, there is a discrete spatial symmetry which 
relates flavors (1,1) and (2,2) to the x and y directions 
in such a way that 



E 2 2 = A'(k) (joko + «AC27l fc l + V F (lj2k2) 



(5) 



D^ u . To proceed we first integrate out the fermions and 
expand the effective action to one-loop order 



ArffK] = (11$ +n ltv )a ll (k)a v (-k) , 
where the bare gauge field stiffness is 
1 
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Now we introduce the dual field = £p V \q v a\, which is 
related to the physical fluctuating vorticity in the theory 
of Ref. 0. We are free to integrate over 6 M with the 
restriction that it is transverse (k^bp — 0). Note that 



£off [bp] = Xobl + Xibj + Xib\ 



(10) 



where {xp} are functions of k^\ 

1 N ^ g^gl x 
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(11) 



At low energies we can neglect the non-divergent bare 
stiffness and thus we set 1 /e 2 = in the above expression. 

The expression ([Io|) is manifestly gauge invariant and 
has the merit of not only being quadratic but also di- 
agonal in the individual components of bp which greatly 
simplifies the computation of the b^ correlation function: 



In the computation of the fermion self-energy, this dis- 
crete symmetry allows us to concentrate on a particular 
pair of nodes without any loss of generality. 

Next, we turn to the gauge field propagator. We first 
work in the Lorentz gauge (fc M a M = 0) and then extend 
our results to a general covariant gauge. To one-loop 
order the fermionic "screening" of the gauge field is given 
by the polarization function 
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^Tr[GUqh;GlS(q + k)^] 
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The above expression can be evaluated by observing that 
it reduces to the isotropic II MI/ (fc) once the integrals are 
properly rescaled (J]. The result is: 
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where we have taken the advantage of the "nodal" met- 
ric (0). This expression is explicitly transverse, i.e. 
kfJlp V {k) = ~n. llv (Jt)k u — and symmetric in its space- 
time indices. It also properly reduces to the isotropic 
expression when vp — i>a = 1- 

As opposed to the isotropic case, it is not quite as 
straightforward to determine the gauge field propagator 



Xp XpXv \,Xi 



(12) 



The repeated indices are not summed over in the above 
expression. 

After this little trick with the integration over bp, it is 
now quite simple to compute the propagator for the orig- 
inal gauge field and in the Lorentz gauge we obtain: 

Dpv{q) = {a^a u ) = e^e v x p -^-(b0bp) . (13) 

By employing the transverse character of {bpb v ) (which 
is independent of the gauge) the above expression can be 
further reduced to 
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(b 2 ) - (bpK) 



(14) 



We use the above result to define a general "covariant" 
gauge for the anisotropic theory as 



Dp,M 



(15) 



where £ is a continuous gauge fixing parameter. This ex- 
pression is justified by the Fadeev-Popov procedure ap- 
plied to the Lagrangian 



Cett[ap] = Ti-p V 



1 2/c kpki-i, 



ap{k)a u {-k) . (16) 
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Note that (b 2 ) can be determined without ever consider- 
ing the gauge fixing terms. The expression ( |l5| ) is our fi- 
nal result for the gauge field propagator in an anisotropic 
"covariant" gauge. 

Having determined the free fermion and screened 
gauge field propagators of the anisotropic theory, we can 
now compute the Dirac fermion self-energy generated by 
the photon exchange to leading order in I/TV: 

j^l^iq-k^D^k) , (17) 

where n is the node index. After some tedious algebra 
this can be manipulated into: 



A 



y/Qc,92pQ0 



(18) 



Here A is the ultraviolet cutoff and the coefficients 77 M are 
functions of the bare anisotropy which can be reduced 
to quadratures. In case of weak anisotropy (vp = 1 + 
6, va = 1) to second order in 5: 
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l-§e-i(4O-70* 



(19) 
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■53j?( 1 -|f + l 4 -S( 4, -'8'"l 



. — (l-7£U 2 ). (21) 

3tt 2 N V 2^ 5 35 V w 1 



In the isotropic limit (vp — v& — 1) we regain 77™ = 
-8(1 - §£)/3tt 2 7V as previously found by others g, ||. 

We are now in position to turn to our main concern: 
the effect of anisotropy at the above critical point of 
isotropic QED3. Before plunging into formal analysis, we 
first make some general physical observations regarding 
the RG flow of the anisotropy. First, by examining the 
Eq. (18) it is clear that if 77" = 7/2 then the anisotropy 
does not flow and remains equal to its bare value. This 
would imply that anisotropy is marginal and the theory 
is described by some anisotropic fixed point. For this 
to happen, however, there would have to be a symme- 
try that preserves the equality 77" = 77J ■ In the isotropic 
QED 3 the symmetry that protects the equality of t^'s 
is the (Euclidean) Lorentz invariance. In our anisotropic 
case this symmetry is broken and thus we generically 
have 77™ ^ 77J . Therefore, the anisotropy runs in the RG 
sense and flows away from its bare value. If we start with 
ap, > 1 and find that 772 1 > 77J 1 at some scale p < A, the 
anisotropy is marginally irrelevant and decreases towards 
unity as we move toward infrared. On the other hand, if 
772 < ?7i j then the anisotropy becomes marginally rele- 
vant, continues increasing beyond its bare value and the 



theory ultimately flows into some new, anisotropic criti- 
cal point. 

The above arguments concerning 772 — 77™ and RG flows 
are on solid ground physically only if they can be made in 
a gauge-independent way. This condition appears com- 
promised by the fact that 77™ 's are gauge dependent quan- 
tities and explicitly include the gauge fixing parameter £. 
However, the difference 772 — 77" is itself gauge-invariant. 
This is seen directly from Eqs. ( p9| - pT| ) where the £ depen- 
dence of all 77's is exactly the same. This cancellation of £ 
dependent terms in 772 — 77™ occurs not only for 8 <C 1 but 
is the general feature of 77™ 's to all orders in anisotropy 
and for any choice of "covariant" gauge fixing. Therefore, 
the RG analysis that follows is fully gauge invariant as it 
should be. 

The renormalized two-point vertex function is related 
to the "bare" vertex via a fermion field rescaling factor 

Zip as = Z^T^ 2 \ It is natural to demand that at 

( 2) 

some renormalization scale p, T R (p) have the form (at 

nodes 1 and 1): T^(p) = 70P0 + 71P1 72P2, where 
Vp and v^ are the renormalized velocities. The above 
equation corresponds to our renormalization condition 
through which we can eliminate the cutoff dependence 
and compute the RG flows. 
To order l/N we can write 



(22) 



where we have used the self-energy ([l8|). Multiplying 
both sides by 70 and taking the trace determines the 
field strength renormalization: 



r » 1 - rS In - 

1 + Vo ln f P 



(23) 



We can now determine the renormalized Fermi and gap 
velocities: 

V l » (1 - In £)(1 + rtf In ^) « 1 - - 7?) k - 
vp p p p 

(24) 



and 

V -l » (1 - 770 11 In A)(i + n f In » 1 - {r,f ^) In A. 

VA P P P 

(25) 

The corresponding renormalized Dirac anisotropy is 
therefore 



old 



i-(% n -??r)in 



(26) 



The RG beta function for the anisotropy is given by: 

ctDiril 1 ~ V?) ■ (27) 
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FIG. 1: The RG /3-function for the Dirac anisotropy in units 
of 8/3ir 2 N. The solid line is the numerical integration while 
the dash-dotted line is the analytical expansion around the 
small anisotropy (see Eq. ( ]19r|2l| )). At qd = 1, (ia D crosses 
zero with positive slope, and therefore at large lengthscales 
the anisotropic QED3 scales to an isotropic theory. 



In the case of weak anisotropy (vp = 1 + 5, va = 1) 
the above expression can be determined analytically as 
an expansion in <5. Using Eqs.(^0[j2l|) we obtain 



Pa 1 



3tt 2 N \ 5 



6(1 + S){2 - S) + 0(S 3 ) 



(28) 



Note that this expression is independent of the gauge 
parameter £. For < S <C 1 the j3 function is positive 
which means that anisotropy decreases in the IR and 
thus the anisotropic QED3 scales to an isotropic QED3. 
For — 1 <C S < the (3 function is negative and in this 
case arj increases towards the fixed point 05 = 1, i.e. 
again towards the isotropic QED 3 . Note that for S > 2, 
P < which may naively indicate that there is a fixed 
point at 6 = 2; this however cannot be trusted as it is 
outside of the range of validity of the power expansion of 
rjfj,. The numerical evaluation of the quadrature for 77^ 
shows that, apart from the isotropic fixed point and the 
unstable fixed point at old — 0, [3 aD does not vanish (see 
Fig. |l|). This indicates that to the leading order in 1/N 
expansion, the theory flows into the isotropic fixed point 
and relativistic invariance is dynamically restored in the 
IR limit. 

Although S„ (q) is a gauge dependent quantity we now 
show that it can still be put to good use in helping us 
extract the physical, gauge-invariant information about 
fermion propagation. To this end, we employ an impor- 
tant result derived by Brown [|l3| in the context of QED4 
which relates G(x — x') = (tp(x)'tp(x')) to the gauge- 
invariant propagator G(x — x'): 



G(x-x') = e- F[x - x '^G(x-x') 



(29) 



where 



and 



F = - I dzdz'J^z)D^{z - z')J v {z') (30) 



J M (z) = (x-x') fl da5 d (z-x' -a(x-x')) . (31) 
Jo 



The exponential factor in (|29j) is just the expectation 
value of the straight line integral of the gauge field 
(exp(z ll ds^a^)) averaged over @. 

The expressions ( p9f]3l| ) hold in an arbitrary covariant 
gauge and reflect the fact that the part of £ e ff[<v] ob- 
tained after integrating out the fermions must be purely 
transverse and thus the longitudinal, gauge dependent 
part (|l6[), enters only at the quadratic order. The lead- 
ing long wavelength dependence is then given by Eqs. 
( f29| - |3~l| ) . Using the long distance scaling of the gauge 
field propagator D flv (s(x — x')) = s 1 ~ d D^ v (x — x') and 
dimensional regularization ]l3| it is straightforward to 
show that F(r) = (d _ 2) 1 (d _ 3) r- A1 £' Att/ (r)r^ , where D^ v (r) 
is the Fourier transform of ( |l5| ) to the real spacetime. 
Computing this expression at the isotropic fixed point in 
the limit d — > 3 we find 



F(r) = 



4(5 



Ntt 2 



„3-d 



3-d 



4(g~2) 
Nn 2 



log(Ar) 



1 



3-d 
(32) 



If £ = 2 (Yennie gauge) F vanishes and the physical 
gauge- invariant propagator G(x — x') coincides with the 
gauge- variant one G(x — x'). An important observation 
follows from the above results: by inserting £ = 2 into our 
expression for the gauge- variant self-energy fllg-21) we 
obtain the anomalous dimension exponent of the physi- 
cal gauge-invariant Dirac fermion propagator in isotropic 
QED 3 , r] = 16/37r 2 iV (jj - 0.27 for N = 2) g ||. 
This exponent serves as the signature of the interact- 
ing infrared critical point which regulates the low-energy 
physics in the chiral symmetric phase ( AFL of Ref . Q] ) 
of both isotropic and anisotropic QED3. 
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